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The effect of Polyakov loop on the QCD phase diagram at high density is studied within 
the Nambu-Jona Lasinio model with Polyakov loop (PNJL model). We point out that the 
color neutrality is missing in the standard PNJL model at finite density. Moreover, we discuss 
how the color-flavor locked (CFL) phase is to be distorted by the inclusion of Polyakov loop. 

§1. Introduction 

The phase diagram of strongly interacting matter has been a subject of theoret- 
ical/experimental work since the foundation of Quantum chromodynamics (QCD). 
The perturbative QCD can be of some help in the extremely high density regime, but 
it is no longer reliable at density of physical interest. The lattice QCD is a powerful 
tool to study such a strongly interacting regimes. However there is a well-known dif- 
ficulty in simulating QCD on lattice at finite density. So exploring phase structure at 
intermediate density, where neither lattice simulations nor perturbative calculations 
can be trusted, remains the subject of various model studies which mimic some of 
basic features in QCD. The Nambu-Jona Lasinio (NJL) model is one of themP and 
it nicely predicts the chiral restoration of QCD at extreme conditions. 21 

The main defect of the NJL model had been the lack of the notion of the con- 
finement. To improve this point, Fukushima included the Polyakov loop dynamics 
into the NJL model,® and the model is now called "Polyakov-Nambu-Jona Lasinio" 
(PNJL) model. This model has two order parameters, qq for the chiral restoration, 
and the Polyakov loop <P for the deconfinement. Even though these two serve as 
the exact order parameters only in the different limits, (m q — > and m q — > do), the 
model enables to interpret nicely some bulk properties of matter observed on the 
lattice on the field theoretical ground.® 

In this work, we will extend the application of PNJL model to color supercon- 
ducting phases at high density 5). In particular, we are interested in: (a) how the 
phase structure in (T, m 2 ///)-plane will be modified by the inclusion of the Polyakov 
loop, and (b) what is the consequence of imposing the color /electrical neutralities 
on the PNJL model with and without diquark condensations. The purpose (a) is 
regarded as the extension of the earlier work [6]), [7]), while (b) is considered as the 
extension of Refs. [8]), [9]). We note there is a parallel development on the role of 
electrical neutrality in the PNJL model at low density.^ 
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§2. The model 

The Wilson line operator is a key quantity whose expectation value plays a role 
of order parameter for deconfinement transition in a pure gauge theory. It can be 
expressed by the background Euclidean temporal gauge field A±(t, x) = igA^r, x)T a 
(T a = 4^; {X a } are the standard Gell-Mann matrices for SU(3) C ) as 

-l/T 



L Q = P T exp [i f lfl drAt(r, x)J . (2-1) 

The Wilson line in the anti-triplet representation can be defined as Lq = LL. In a 
pure gauge theory with zero chemical potential for quarks (fx = 0), these are regarded 
as the operators associated with heavy (anti)quark excitation in the gluonic heat 
bath at temperature T; loop <P = -^-(trLg)T (anti-triplet loop = -^(trLg)y) is 
related to the Free energy of single (anti)quark excitation in the gluon medium by 
<P = e~ FQ l T and ^ = e~ F Q^ T '. In this case, one finds eventually = <P = #*. In 
the full gauge theory with dynamical quarks and with a finite chemical potential for 
quarks, however, this is no longer the case because quark and antiquark propagate 
differently in each direction of imaginary time. The detailed analysis of a matrix 
model shows that ^ and <P certainly differ from each other but still both stay real; 
this is due to the imaginary piece in the action proportional to the imaginary part of 
■j^-tiljQ which comes from integration of dynamical quarks and is C-odd quantity.^ 
In the PNJL model, the loop ^ can be parametrized by eight real parameters 
{<pi, (f2, • • • , fs} each of which has a dimension of energy, as 

<2>L4 4 ] = ^tre^/ T , M = EaU 1 P«T a . (2-2) 
Using this A^ field, the PNJL model is given by the following Lagrangian density 
C[q, q; A 4 ] = q{i{f[M] + 7o(a* + SfX eS ))q + %qP vv >q T q T P m >q ~ U{T, <P[A 4 }). (2-3) 

q = (qur,qd g ,q S b, q ug ,qdr, q S r,q U b, qdb,q sg ) T is the quark field. £> M = <9 M - S^A^ is 
the covariant derivative through which the Polyakov loop can change the nature of 
propagation of dynamical quarks. G parametrizes the strength of attractive coupling 
in the color-flavor channel P m > = C^e^ijt^ab {Prrrf = 70-P,L/7q)- In this work, we 
take the CFL type diagonal ansatz for diquark condensation,^ i.e., 




{q T P w q} = A 2 =A W . (2-4) 



r\ (rf) stands for the flavor (color) index. We work within the chiral SU(2) limit set- 
ting m u = m<f = 0, and take into account the strange quark mass m s within the high 

2 

density approximation. So we set <5/x e ff = — l^eQ + ^T^+ ^T^, — ^-diag.(0, 0, l)f x l c , 

where Q = diag.(2/3, -1/3, -l/3) f x l c , T 3 = l f x |A 3 , and T 8 = l f x -i-Ag. U(T,$) 

V 3 

is the Polyakov loop potential which controls the confinement /deconfinement tran- 
sition in the pure gauge sector, whose detailed form will be given later. 
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The real part of effective potential within the high density effective theory 
(HDET) comes out to be 

Mf2(A v ,<p a ) = U(T,$)-^-^-2g± 

-ELl f^^ + /^&2Tln(||l + e-^^)/ r || " 



12tt 2 1 J 2vr 2 
^ V^9 f^c (/•»+i||) 2 <«|| 



,A 



(2-5) 

2 1| is the quark momentum measured from the Fermi surface p = \i. uj c is the ultra- 
violet cutoff needed to regularize the divergent third line. The complex energies 
{Ei, Ei,--- , E§\ are defined by choosing the eigenvalues of non-hermitian matrix 

V A^aberjij — Z|| + <5^* ff — %A\ J ' 

such that Ea — ► — M^gl when — > is satisfied. We use the cutoff (ui c ) depen- 
dent coupling ^ = 2t^_ In ( 2 i/"^ o ) where Z\ is the magnitude of CFL gap parameter 

at T = in the chiral limit (m s = 0). With this convention, the effective poten- 
tial is only weakly (logarithmically) divergent, so the gap equations and neutrality 
conditions have well-defined finite limit as lo c — > oo. 

Parameter reduction via gauge invariance: The effective potential is a function 
of three gap parameters {A^} and eight parameters {if a } f° r the Wilson line matrix. 
In certain cases, the gauge invariance is helpful to reduce the dynamical variables. 
A„ = is such a case; the number of parameters for can be reduced from N% — 1 
down to N c — 1 as we see below. Since the Wilson line transforms as Lq — > gLng~ l 
by the gauge transformation with g being the arbitrary SU(3) matrix, we see 

Q(A A ) = Q(gA 4 g- 1 ), (2-7) 

from the gauge invariance. We can always make A A diagonal by choosing the suitable 
g as gA^g -1 = foTs+cfrsTs, so we can work with this simplified ansatz for A\ without 
loss of generality. The effective potential as well as the Polyakov loop ^ is a function 
of {4>3, 4>s} i n this case. This procedure is just like the change of integration variable 
from SU(3) to its eigenvalues by integrating out six phase variables. However, once 
diquarks come into the problem, this simple reduction does no longer work. In fact, if 
we try to diagonalize A4, the diquark condensate also suffers from the gauge rotation 

A A -> gA A g-\ A m> -» A' w = (Ag- 1 ) . (2-8) 

We note that the new condensate matrix A 1 , is no longer restricted to be of di- 
agonal form in the color-flavor space. A A can be diagonalized by g which can be 
parametrized by six "phase" parameters; at the same time the diquark condensate 



*' We have also an imaginary part of the effective potential, which may viewed as a sign problem, 
ft can cause a splitting of 4> and <P, but here we simply discard itP Accordingly, we have $ = <P. 
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acquires this phase rotation, so the new diquark condensate is to be parametrized 
by nine parameters. Thus in principle, we have two choices; one is (1) to work in the 
diagonal form of A4 with the generalized off-diagonal ansatz for diquark condensate. 
The other is (2) to work in the standard diquark ansatz with A„ but with a general 
A4 parametrized by eight parameters. These proper treatments require us to deal 
with all eleven variational parameters. Leaving these proper arguments to our future 
plan, we here work within the simplified ansatz for the ground state that the diquark 
condensate is diagonal {A^} even after fixing the gauge which diagonalizes A4; this 
means we take only the diagonal entries of A4, {^3, tps}, as the variational variables. 

Although the continuous gauge freedom should be considered to have gone away 
to diagonalize A^, there remain six discrete gauge transformations each of which 
leaves <P unchanged; these are the elements of permutation of fundamental color 
indices. Also in accordance with discarding 9/2, we further put ^ to be real. For 
this, we take ips = 0, so in this way, the gauge is completely fixed in our calculations. 

Finally, we specify the Polyakov loop potential. We adopt the following formUj) 

U^m = _^)0*0 + ft ( T ) log (l _ g 0*0 + 4(^*3 + 03) _ 3(0*0)2) ; (2-9) 

with b 2 {T) = ao + ai (^) +a 2 (^) 2 , b(T) = b 3 (^) 3 . T is the value of the transition 
temperature for deconfinement in pure gauge, i.e., To = 270 MeV. The logarithmic 
term was first proposed in .3) and it is nothing but the Vandermonde determinant, 
i.e., the Jacobian associated with the change of dynamical variables from the SU(3) 
matrix to its eigenvalues [T3]) . For more details for parameter setting used in our 
numerical analysis, the readers are referred to Ref. [5]). 

§3. Results and Discussion 

Color density associated with color symmetry breaking: Before going into 
the full calculation, we study the simplified case with unpaired matter in the chiral 
limit, A v = m s = 0, just to illustrate the importance of imposing color neutrality 
in the PNJL model at finite fi. In Fig. [H we show the color T% density, (q^Tgq), 
and the Polyakov loop $ as a function of T. Surprisingly, Ts color density takes 
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Fig. 2. (Reprinted figure with permission from [H. Abuki et al, Phys. Rev. D 77, 074018 (2008)] Copy- 
right (2008) by the American Physical Society), (a): Phase diagram in (-Sj , T)-plane at Ao = 
60MeV, fj. = 500 MeV. (b): The same as (a) but without the Polyakov loop. 

nonzero valued It can be shown that this is the case except for two different limits, 
T —* and T — > oo (<P — > 1).® The Polyakov loop <P is the colorless object, so one 
might think it is strange to have nonvanishing color density. The reason is simple; 
we are breaking color symmetry in addition to Z% center symmetry by introducing 
the constant A4 background field. In fact, <P = -^trLg is invariant under the color 
rotation Lq — > gLqg^ 1 ; it changes its value only under color transformation which 
is not exactly periodic in imaginary time but only up to Z3, i.e., <P — > z<P under 

Lq -> gLQ(zg)- 1 , where z e Z 3 . (3-1) 

Since we assumed the constant background to parametrize <P, and A± in contrast 
to <P is not invariant under color rotation, we have broken the color symmetry in 
addition to Z3 symmetry. This is unexpected, undesirable feature of the PNJL model 
and may be considered as the model artifact. It could be dangerous for theoretical 
foundation of the model itself, but we do not discuss further this problem here. 
Instead, we simply assume that the model is still useful once we impose the vani shing 
color density as the constraint by tuning color chemical potentials {^3, fJ-s} — 
The phase diagrams: The phase diagram coming out from our high density PNJL 
model is displayed in Fig. [2] (a). For comparison, we have also shown in Fig. [2j (b), 
the phase diagram calculated with the model without the Polyakov loop. From these 
figures, the impact of the Polyakov loop dynamics on the quark Cooper pairing is 
clear; it has two major effects. 

(a) First, we notice that the Polyakov loop dynamics stabilizes the 2SC phase 
significantly. In fact, the critical temperature for the 2SC-to-unpaired phase tran- 

*' It should be noted that the color density itself is a gauge dependent quantity and thus should 
depend on the choice of the gauge. With our diagonal representation of A4, with ips = 0, the T% 
color density becomes finite as we observed above. If we selected a different gauge, the other entries 
of octet color density {(q^Taq)} should have appeared. The important thing is, however, whichever 
gauge we choose, some color density should become finite; in fact the squared sum of the octet color 
densities is shown to be the gauge independent quantity 1 14|) . 

**' We checked that off-diagonal color charge densities (q^T a q) (for o / 3, 8) are automatically 
vanishing for all the situations we are interested in here. 
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sition at m s = is almost doubled by inclusion of the Polyakov-loop dynamics. 
This point can be understood by the observation that the Polyakov loop suppresses 
the thermal excitation of colored quasiquarks which tend to break the Cooper pair 
condensate. 5 - Numerically, the factor of enhancement of T c is 1.8 which is in good 
agreement with our analytical estimate 1.79.® 

(b) Second, as a consequence of the effect (a), we have the color-flavor unlocking 
transition even at m s = 0. One may wonder why the SU(3) flavor symmetry should 
be broken down to the isospin SU(2) in the (u, <i)-sector, and why not either in 
(s, u) or (d, s) sector. This is strange because at m s = the flavor SU(3) symmetry 
is perfect so how can the flavor-blind Wilson line distinguish them? Actually, the 
fact that we have the isospin symmetry intact is directly attributed to our model 
assumption mentioned in Sec. [2j we are limiting ourselves to treat only two out of 
eight parameters for the SU(3) C matrix, Lq. As noted, in principle this can not be 
justified in our case because the gauge is already fixed in the diquark sector once we 
put the ansatz for diquarks to the diagonal form A^i = diag.(Z\i, A2, ^3); therefore 
there remains no continuous gauge freedom to rotate Lq to a diagonal form. So 
this limitation should be rather viewed as an ansatz for the many possible ground 
states, such as the color-flavor locking ansatz. For the proper treatment, we should 
take into account all the eight parameters to represent the Wilson line Lq. It is 
possible after making such a proper treatment, either that the ground state prefers 
the diagonal form of Lq or that the ground state we obtained here turns out to be 
one of several degenerated ground states. We defer this task in future. 
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